Diverse corrugation pattern in radially shrinking carbon nanotubes 
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Stable cross-sections of multi-walled carbon nanotubes subjected to electron-beam irradiation are 
investigated in the realm of the continuum mechanics approximation. The self-healing nature of sp 2 
graphitic sheets implies that selective irradiation of the outermost walls causes their radial shrinkage 
with the remaining inner walls undamaged. The shrinking walls exert high pressure on the interior 
part of nanotubes, yielding a wide variety of radial corrugation patterns (i.e., circumferentially wrin- 
kling structures) in the cross section. All corrugation patterns can be classified into two deformation 
phases for which the corrugation amplitudes of the innermost wall differ significantly. 



PACS numbers: 61.46.Fg, 46.70.De, 62.50.-p, 61.80.Az 
I. INTRODUCTION 

Carbon nanotubes exhibit remarkable flexibility when 
subjected to cross-sectional deformation^— Such flexi- 
bility has been evidenced by spectroscopy and diffrac- 
tion measurements in which hydrostatic pressure on the 
order of a few GPa caused flattening and polygoniza- 
tion in the cross section4r— Molecular-dynamics (MD)2r— 
and density-functional-theory-based simulation s 10 i 11 re- 
vealed more details about the mechanism of radial defor- 
mation under pressure. Radial deformation can also be 
caused by a nanoindcntation^ under which the elastic 
response of single-walled nanotubes (SWNTs) is deter- 
mined by a universal law that depends on the tube radius 
r and applied load^ As the radial deformation of nan- 
otubes is strongly correlated with their electroni c 12 ' 14 " — 
and optical properties, its clarification should be im- 
portant from the viewpoint of the use of nanotubes as 
components of nanoscale devices. 

The cross section of an isolated SWNT exhibits a 
circular-to-clliptic transition at a critical hydrostatic 
pressure p c . The value of p c decreases with r a o 18 i 19 
p c oc r -3 , which agrees with the conclusion of the clas- 
sical continuum theory^ This scenario, however, fails in 
the case of multi- walled nanotubes (MWNTs). It has 
been shown that MWNTs consisting of several tens of 
concentric walls undergo a novel cross-sectional deforma- 
tion, called radial corrugation,— in which the outer walls 
exhibit wavy structures along the circumferential direc- 
tion. The radial corrugation originates from the multi- 



layered nature, i.e., the competing effects between the 
mechanical instability of outer walls with large radii and 
the radial rigidity of inner walls with small radii. 

The corrugation behavior of a MWNT is believed to 
change if it occurs in an embedding elastic medium. 
Suppose that a MWNT is embedded in a large elastic 
medium that is in complete contact with its outermost 
wall. A uniform shrinkage of the medium exerts an ex- 
ternal force on the outermost wall, as similar to the case 
of hydrostatic pressure. But a difference arises when the 
medium possesses moderate stiffness, since the energy 
required to deform the medium should be responsible 
for determining the stable corrugation pattern. In fact, 
when a thin flat plate is embedded in an elastic medium 
and comressed, it shows a uni-directional corrugation 
whose profile is strongly dependent on the medium's 
stiffness! 22 ' 23 This fact implies an uncovered corrugation 
mechanism peculiar to embedded MWNTs, it i.e., din- 
stinct from the hydrostatic pressure cases, which imposes 
new challenges. 

In this article, we demonstrate that the presence of a 
surrounding elastic medium triggers diverse variations of 
corrugation modes that cannot be observed in MWNTs 
under hydrostatic pressure. These diverse modes are 
found to be grouped into two corrugation phases that 
exhibit a significant difference in the corrugation am- 
plitude of the innermost wall of the MWNT. These 
findings will help to tune the core-tube geometry of 
MWNTs, thus providing useful information for devel- 
oping nanofluidio^ir— or nanoelectrochemical^l devices 



whose performance depends on the geometry of the in- 
ner hollow cavity of nanotubes. 



II. METHODOLOGY 
A. Self-healing nature of graphite sheets 

The current work has been inspired by the discovery 
of the self-healing properties of carbon nanostructures 
under electron-beam irradiation.— ~— Banhart et a/.— 
found that carbon onions, assemblies of concentric spher- 
ical carbon shells, shrink toward smaller onion struc- 
tures when irradiated. This shrinkage is attributed to 
the knock-on collision of carbon atoms from constituent 
shells; the resulting vacancies are healed by annealing 
reconstruction 3 - 3 - that shrinks the system with reduced 
intershcll spacings. As a result, extreme pressure arises 
in the core of spherical carbon shells. Using a similar irra- 
diation technique, Sun et al£^ synthesized MWNTs with 
reduced interwall spacing and revealed high pressures of 
several tens GPa within the hollow cavity. It is thus ex- 
pected that self-contraction of only outermost walls of a 
MWNT (possibly realized by finely-focused beam) will 
provide the setting for radial corrugation of the interior 
part of the MWNT. 

It should be borne in mind that in the carbon onion ex- 
periments, diamond nucleation occurs at the center of the 
onion as the induced pressure is too high for the original 
sp 2 bonds to persist To make the current work con- 
vincing, therefore, it is important to examine whether 
the critical pressure p c for the corrugation is low enough 
to prevent the breaking of sp 2 bonds within the inner 
walls. We will revisit this issue in Section III B. 



B. Continuum model of radially shrinking MWNT 

When treating many-walled nanotubes, atomistic sim- 
ulations are realistic and accurate but they demand huge 
computational cost in general. Thus we used a simplified 
model based on the continuum approximation! 18 i 34 ~ ^ 
First, we map an MWNT onto concentric cylindrical 
shells endowed with^i the in-plane stiffness C = 345 
nN/nm, the flexural rigidity D = 0.238 nN-nm, and Pois- 
son's ratio v = 0.149, in which the inter-shell separation 
is defined by 0.3415 nm. Figures HJa)-(c) illustrate the 
shrinking process considered in this paper. We assume an 
(AI + A)-wallcd nanotube under irradiation that kicks off 
a portion of carbon atoms located within M N) outer 
walls surrounding N inner walls. Then, self-contraction 
of the M walls results in high pressure that exerts on the 
N walls that they encapsulate. For the sake of analyti- 
cal arguments, we simplify the M walls by a continuum 
elastic medium and assume the outermost boundary of 
the 7V-walled tube to be in contact with the inner sur- 
face of the medium throughout the shrinking process (see 
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FIG. 1: (color online) Illustration of radial shrinking hypoth- 
esis of MWNTs. (a) An (M + A)-walled nanotube is sub- 
jected to electron-beam irradiation, (b) Irradiation with fine 
resolution kicks off a portion of carbon atoms located within 
A/(2> N) outermost walls, (c) The eroded M walls are healed 
and radially contracting, as a result of which extreme pressure 
is exerted on the N innermost walls, (d) Typical corrugation 
modes of the pressed A-walled nanotubes are shown with the 
mode index n. 



Figs, [ljb) and (c)). As a consequence, various corruga- 
tion patterns arise in the cross section of the A-wallcd 
nanotube. Only a few examples of the corrugation pat- 
terns are shown in Fig. [TJd) together with the associated 
mode index n representing the wave number in the cir- 
cumferential direction. 

It is noteworthy that the mechanica l 38 ' 39 and 
structural^ consequences of irradiation in few-walled 
nanotubes have already been explored by atomistic 
simulations. Xu et at numerically reproduced the 
irradiation-induced high pressure within a double- walled 
nanotube 3 ^; they artificially removed carbon atoms from 
the outer wall and observed the subsequent self-healing 
process that causes an extreme pressure acting on the 
inner wall. A similar healing process is assumed to oc- 
cur when the outer walls of a many-walled nanotube 
are eroded locally, as shown in Fig. QJb) . It should be 
noted that in the existing experiment, it is possible to 
remove carbon atoms from an MWNT with monolayer 
precision^ This fact supports the validity of our hy- 
pothesis that a part of the carbon atoms within only a 
desired number of outermost walls is removed locally. 



C. Energy formulation 

The stable cross-sectional shape of the embedded tube 
is obtained by minimizing its mechanical energy U per 
unit axial length^ 

u = u D + Ui + u M + n. (1) 
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The first term Ur> = J2iLi ^d' with the definition 
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represents the deformation energy; and Ki are, respec- 
tively, in-plane and bending-induced strains of the ith 
wall having the radius n and 9 is a circumferential coor- 



by the decomposition Ui(p, &) = uf\p) + 6ui{9), where 

u\ (p) oc p describes a uniform radial contraction at 
p < p c and 5ui(9) describes a deformed (non-circular) 
cross section observed immediately above p c . Similarly, 
we can write Vi(p,9) — Svi(9), because u|° '(p) = at 
p < p c . Applying the variation method to U followed by 
the Fourier series expansions 



dinate. The second term Ui 
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van der Waals (vdW) interaction energy of all adjacent 
pairs of walls. Here, Ui is the radial displacement of the 
ith wall, and the coefficients c^- are derived through a 
harmonic approximation of the vdW forced accociated 
with the vdW potential V(r) = 4a[(/3 / 'r) 12 - (P/r) 6 } 
withH a = 2.39 meV and ,5=0.3415 nm. The final term 
in Eq. ([!} is the negative of the work done by p dur- 
ing cross-sectional deformation. It can be proved that^i 
all the three terms are functions of Ui(p,9) and the cir- 
cumferential displacement Vi(p,9) of the ith wall under 
P- 

The remaining term Um in Eq. ([1} is the elastic en- 
ergy of the surrounding medium. To derive it, we as- 
sume that the medium is homogeneous and isotropic with 
Young's modulus Em = 100 GPa (that corresponds to 
the modulus of amorphous carbo n 44 ' 45 ) and Poisson's ra- 
tio um = v. In polar coordinates, the radial and circum- 
ferential components of normal stress in the medium are 
denoted by oy and ag, respectively, and the shear stress 
is denoted by r r g; all the three quantities are functions 
of r, 9. Then, Um is determined by oy and T r g at r = 
as 



l) 2 d9 accounts for the 8uj(8)=y] 8pg{n) cos n9, Svj(9)=y^ j 5vj{n) sinnfl, (6) 



U M = U$ + AU$ ) 



u 



(o) 

M 



AU 



(n) 
M 



2 jo 
2 ,/i 



r(°) 



-(») 



u [ ^d9, 



Sun 



» 



(3) 
(4) 

5v N ) d6, (5) 



where Sun and 5v n describe the corrugation amplitudes 
of the outermost wall of the embedded MWNT; see 
Eq. ©. The superscripts (0) and (n) indicate that the 
quantities correspond to a uniform contraction and radial 
corrugation, respectively. In plane words, llff repre- 
sents the energy required for uniform radial contraction 
of the MWNT remaining in contact with the medium, 
and AUm does for radial corrugation with the mode in- 
dex n. Details of the derivation of Um are presented in 
Appendix A. 



D. Corrugation mode analysis 



we obtain the matrix equation^! Mu = 0; the vector 
u consists of Sp,i(n) and 5vi(n) with all possible i and 
n and the matrix M involves one variable p and other 
material parameters. Finally, we solve the secular equa- 
tion det(TW) = with respect to p to obtain a sequence 
of discrete values of p among which the minimum one 
serves as the critical pressure p c . Immediately above p c , 
the circular cross section of MWNTs becomes radially 
deformed, as described by 



u.i(9) = uf\p c ) + Sfii(n) cosn9, 
Vi(9) — 8vi(n)sa\n9, 



(7) 

(8) 



where the value of n is uniquely determined by the one- 
to-one relationship between n and p c . 

We restrict our attention to the elastic limit (i.e., dfii 
and 5vi be infinitesimally small); hence, neither delam- 
ination nor strong stacking effects are considered in the 
subsequent arguments, although these effects may change 
the corrugation patterns.. 



III. RESULTS 
A. Cross-sectional view 

Figure [TJd) illustrates a cross-sectional view of typical 
deformation modes: the two left-hand-side panels show 
"inward-deformation" modes with radial-corrugation 
mode indexes n = 4 and 5, and the right-hand-side 
panel shows an "outward-deformation" mode with n = 
9. In the inward-deformation mode, the innermost 
walls exhibit significant corrugation amplitudes as com- 
pared to the outside walls. Conversely, in the outward- 
deformation mode, the innermost wall maintains its cir- 
cular shape. Which class of modes is observed immedi- 
ately above p c depends on the values of the innermost 
tube radius R(= and N under consideration. As 
shown below, larger R and smaller N favor the inward 
mode with larger n. 



Our objectives are to determine: (i) the optimal dis- 
placements Ui and Vi that minimize U under a given 
value of p, and (ii) the critical pressure p c above which 
the circular cross section of an MWNT is elastically de- 
formed into a non-circular one. These are accomplished 



B. Critical pressure for radial corrugations 

Figure [5] (a) shows the TV-dependence of p c for various 
conditions of R. For all Rs, p c exhibits two shallow peaks 
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radii grow with N, implying the occurrence of outward- 
deformation modes at large N. 

It also follows from Fig. [U^a) that a large portion of 
PcS data lies on the order of several GPa, though they 
occasionally exceed ten GPa or more at certain limited 
conditions. In the latter conditions, the continuum shell 
approximation may break down due to interwall sp 3 /sp 2 
hybridization bonds generated by pressure. In fact, ear- 
lier MD simulations^ suggested that radial pressure of 
approximately 40 GPa (or 30 GPa) destroys the tubular 
structure of an encapsulated SWNT having the radius 
~ 0.7 nm (~ 1.2 nm) through the formation of hybrid 
sp 2 /sp 3 self linkages inside the tube. It is thus con- 
jectured that similar hybrid bonds between neighboring 
walls in MWNTs for R > 1.0 nm occurs above a threshold 
pressure comparable to (or maybe less than) tens GPa. 
Precise determination of the threshold requires atom- 
istic calculations; hence we proceed arguments bearing 
in mind that the obtained results are physically relevant 
only when p c is lower than the threshold (more or less 
tens GPa) that remains to be determined. 

Figure shows the three-dimensional plot p c on 

the R-N plane. The plot exhibits a ridge line that ex- 
tends from the top at (R,N) = (1.0,2) to the skirt at 
(R, N) = (4.0, 10) of the p c -surface. The ridge line corre- 
sponds to a phase boundary that separates the inward- 
and outward-deformation phases. In fact, an inward- 
deformation mode occurs in the region to the left of the 
ridge line, as we will find in the phase diagram in Fig. [3] 
(a). 



C. Phase diagram 



FIG. 2: (color online) (a) Critical pressure curve p c (N) as 
a function of the total number of concentric walls N con- 
tained in a MWNT. The innermost tube radius R ranges 
from 1.0 nm to 6.0 nm as indicated, (b) Three-dimensional 
plot of p c (R,N) on the R-N plane. A ridge line extending 
from the top to the skirt of the p c -surface, depicted by a 
bright dashed curve, corresponds to the phase boundary that 
separates the inward-deformation phase from the outward- 
deformation phase (see text for definitions of the two phases). 



(one upward and one downward) whose positions shift to 
larger N with an increase in R. One might expect that a 
larger N leads to a larger p c , because an increase in the 
concentric walls would enhance the radial rigidity. This 
conjecture holds for MWNTs with intermediate values 
of N between the two shallow peaks; for instance, p c 
for R = 4.0 nm increases slowly with N between N 
4 and 10. However, to the right of the upward peak, 
p c decreases monotonically with N, in contrast to the 
conjecture above. Such decay of p c at large N arises 
from the mechanical instability of outside walls whose 



Figure [3£a) shows a phase diagram of the radial corru- 
gation modes in MWNTs observed above p c . A stepwise 
bright line represents the phase boundary between the 
inward-deformation phase (below the line) and the out- 
ward phase (above). Figure [3fb) shows a contour plot 
of p c (R,N) that evidences a strong correlation between 
the ridge line of p c (R,N) (i.e., the solid slanted line in 
Fig. EJb)) and the phase boundary shown in Fig. [3ta). 

A salient feature of Fig. Eta) is the absence of an ellip- 
tic deformation phase (n = 2) within the ranges of R and 
N we have considered. The absence of the n = 2 mode 
is in contrast with the cases of MWNTs and SWNTs 
under hydrostatic pressure; in fact, the n = 2 mode in 
the latter two cases is the primary mode observed in a 
large domain of the R-N space. Furthermore, we have 
obtained a wide variety of corrugation modes for various 
values of R and N, where the variation of n is systematic 
with the changes in R and N. In the inward-deformation 
phase, for instance, larger R and smaller N favor corru- 
gation modes with large n. Contrariwise, in the outward 
phase, smaller R and larger N favor corrugation modes 
with large n. 
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FIG. 3: (color online) (a) Phase diagram of the radial corrugation for embedded MWNTs. A slanted solid line represents 
the phase boundary between the inward-deformation phase (below the line) and the outward phase (above), (b) Contour plot 
of the critical pressure p c (R,N) on the R-N plane. A solid line represents the ridge line appearing in the three-dimensional 
Pc-surface (see Fig. H£b)). 



D. Corrugation amplitudes 

The contrasting difference in corrugation amplitude 
distributions between the inward and outward phases is 
quantified by plotting the deformation amplitudes Sfii 
introduced in Eq. ([7]) . Figure 2] shows the normalized 
deformation amplitudes, £j = | (5/2^/(5/2^1, of individual 
concentric walls for iV-wallcd nanotubes with different 
N and R = 6.0 nm be fixed, ft follows from the figure 
that the value of £i suddenly drops off from £i > 1.0 to 
£i < 1.0 across the phase boundary (i.e., from A^ = 8 
to 9 with R = 6.0 nm being fixed). To understand this 
visually, we show in Fig. [5] the sequential variation of the 
corrugation amplitude distribution in the cross section 
near the phase boundary, In the inward mode (N < 8), 
the restoring force exerted by the surrounding medium 
is relatively strong so that it tends to prevent radial de- 
formation of outer walls; this is why inner walls exhibit 
significant corrugation amplitudes to lower the energy U 
of the system. By increasing TV, such restoring force ef- 
fects become ineffective to support the radial instability 
of outer walls, as a result of which the system falls into 
the outward phase for A^ > 9. A similar scenario applies 
if we fix N and modulate R in the vicinity of the phase 
boundary. 



IV. DISCUSSIONS 

The present results are based on the approximation 
that the radially shrinking part of an MWNT under 
electron-beam irradiation is mapped onto a continuum 
elastic medium with homogeneous and isotropic elastic- 
ity of the modulus Em = 100 GPa. The validity of the 
homogeneity and isotropy hypothesis depends on the fol- 



lowing two effects of irradiation on the mechanical stiff- 
ness of MWNTs: irradiation reduces the axial stiffness 
because it creates vacancieS ) 46 i 47 and simultaneously, it 
enhances the radial stiffness owing to the production of 
covalent bonding between adjacent walls^ A quantita- 
tive examination of the degree to which this hypothe- 
sis holds requires elaborated measurements or large-scale 
atomistic simulations, and therefore, we have no avail- 
able data for the proof. Instead, it is possible to gener- 
alize the theoretical method by considering the possible 
anisotropic elasticity in the medium. Our preliminary 
calculations showed that moderate anisotropy causes lit- 
tle modification in the phase diagram. More detailed 
results will be shown elsewhere. 

We have also performed corrugation analyses by im- 
posing other values of Em than Em = 100 GPa. It was 
found that larger (smaller) values of Em result in more 
(less) number of corrugation modes observed in the phase 
diagram with fixed ranges of R and A?" set in Fig. 
With decreasing Em, the phase boundary shifts down- 
ward until Em = 0.1 GPa, below which the boundary 
disappears leaving the outward phase as the only possi- 
bility. Therefore, our predictions of diverse corrugation 
patterns and the two corrugation phases are possible for 
any choice of Em as far as Em > 0.1 GPa. 

From an engineering perspective, the selectivity of the 
innermost wall geometry by tuning the material param- 
eters R, N and p may be useful in developing nanotube- 
based nanonuidio^— or nanoelectrochemical devices^ 7 -, 
because both utilize the hollow cavity within the in- 
nermost tube. A very interesting issue from the aca- 
demic viewpoint is the effect of the core tube defor- 
mation on the physical and chemical properties of in- 
tercalated molecules confined in the hollow cavity. It 
has thus far been known that various types of interca- 
lated molecules (diatomic gas, water, organic, and tran- 
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FIG. 4: (color online) Normalized corrugation amplitude f j = 
\5pd/5p,N\ of each i-th concentric wall of N- walled nanotube 
with R = 6.0 nm. Inward- (outward-) deformation modes ar 
characterized by the value of £1 larger (smaller) than 1.0. 

sition metal molecules, etc.) can fill the innermost holloa 
cavities of nanotubes^ and exhibit intriguing behavior 
that are distinct from those of the corresponding bul 
systems! 48 ' 49 These distinct behaviors originate from th 
fact that the intermolccular spacings become compara 
ble to the linear dimension of the nanoscale confining 
space. Therefore, the core tube deformation that breaks 
the cylindrical symmetry of the initial nanoscale com- 
partment will engender unique properties of intercalated 
molecules that are peculiar to the constrained condition 
in a radially corrugated space. 

Another interesting implication of our results is a 
pressure-driven change in the quantum transport proper- 
ties of 7r-electrons moving along the radially corrugated 
nanotube. It has been known that mobile electrons whose 
motion is confined to a two-dimensional curved thin layer 
behave differently from those on a conventional flat plane, 
because the geometric curvature of the layer effectively 
yields an electromagnetic field££~— that can affect low- 
energy excitations of the electrons. A quantitative exam- 
ination along with sophisticated atomistic simulations^ 
should reveal novel MWNT applications based on radial 
corrugation. 

V. CONCLUSION 

We have theoretically shown the presence of diverse 
radial corrugation modes in the cross section of MWNTs 
that shrink radially under irradiation. Using the contin- 
uum elastic approximation, we have established a phase 
diagram that enables a desired corrugation pattern to be 
obtained by tuning the innermost wall radius R and the 
total number of concentric walls N. We have also found 
that all corrugation patterns are classified into two de- 
formation phases between which there exists a significant 
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FIG. 5: (color online) Cross-sectional view of the inward- 
(N = 6, 8) and outward- (TV = 9) corrugation modes. 



difference in the corrugation amplitude of the innermost 
wall of the nanotube. We believe that the results provide 
useful information for developing carbon-nanotube-bascd 
devices that utilize the nanoscale hollow cavity within the 
core of concentric carbon walls. 
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Appendix A: Derivation of Um 

In this Appendix, we derive the deformation energy 
Um of an elastic medium surrounding a MWNT. The 
mechanics of an clastic medium is governed by the 
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stress function <f> that satisfies the so-called compatibility 
equation^ 

{dl + r^dr + r^def ${r,6) = 0, (Af) 

where d r = d/dr and do = d/89. Once <j> is obtained, we 
can deduce the radial and circumferential components of 
the normal stress, o r and og, respectively, and the shear 

I 



stress r r g as follows. 

o r = {r^dr + r- 2 d 2 e ) 0, o e = d 2 r <j>, (A2) 
Tr6 = d r (r~ 1 dg)cf>. (A3) 

By definition, the strain components e r ,eg,^/ r e are given 
by the matrix form 
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where u = u(r,&) and v = v(r,9) are respectively the 
radial and circumferential displacements of a volume el- 
ement in the host medium. 

The general solution of Eq. I|A1|) is given by 
4>{r-,0) = Y^=o4>n{r,0), where <f> n {r,9) = f n (r) cos n9 + 
g n (r) sin n9 and^ 

fa = a a logr + b a r 2 logr + c a r 2 +d a + a a + P Q r 2 6(A5) 
fi = a 1 r~ 1 +bir\ogr+ar 3 + dir+ai9, (A6) 
fn = a n r- n +b n r 2 - n +c n r 2+n +d n r n , (n > 2) (A7) 

with similar definitions of g n . The zeroth component cj>o 
represents a uniform contraction of the circular cross sec- 
tion, thus corresponding to the energy 

U$ that we have 

introduced in Eq. (QJ. The first one 4>\ implies a rigid 
body translation that is irrelevant to our consideration. 
Other components 4> n for n > 2 describe radial corruga- 
tions with mode index n, thus providing the energy AUm 
given by Eq. ([5]) . In the following, we set fa = ciq log r 
and /„ = a„r~ n + b n r 2 ~ n instead of Eqs. (|A3|) and (|A7J), 
respectively, in order to obtain physically relevant solu- 
tions of o r ,og, t t b that decay with increasing rJ& 

We now evaluate the explicit forms of and AU^ . 
To derive U^,\ we consider a specific solution of (|A1[) 
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that has the form of 



bo, and then substitute it 



in Eqs. (|A2|) and (|A3|) to obtain a r 
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0. Hence, it follows from Eq. (|A4 



that u (°)(r) = (1 + v M )ao/(E u r), w (0) = 0. Com- 
plete contact between the medium and the outermost 



wall implies «^(rjv) = , and the elastic nature of 
the medium implies (Tr°^|r=rjv = K o u 



(0) 
N 



with the stiff- 



and 



t (0) 



Next, we consider the energy Auffl (n > 2) that cor- 
responds to the radial corrugation of the nth order. We 
only consider the cosine term in 4>n without loss of gen- 
erality, which is based on our assumption of cosine radial 
displacement 8ui{9) (see Eq. (J5|)). A similar procedure 
to the case of n = yields 



,(«). 



■ {—n(n+l)a n r 2 -(n-l)(?i + 2)6„}r "cosn6>, 



(A8) 



■ {n(n + l)a n r 2 + (n- l)(n— 2)6„}r "cosn#, 



(A9) 

r W = {_ n ( n + l)a n r~ 2 -n{n-l)b n }r- n smn9, (A10) 
leading to the results 
+ 



B w (r, 



w (n) (r, 



na„r 



+ {n + 2(l - 2i/M)}M rl ~" C0S "- 6 ' 5 ( A11 ) 

(! + ^m) r _ 2 
= — = \na n r 

+ {n - 4(1 - i/ M )} b„] r 1 "" sinnfl. (A12) 



Parallel discussions to the n = case together with the 
formula 



» 



{knSpLN + /ci2<5i / 7v) cos n(9, (A13) 
(k2i6p.N + ^SPn) sinn#, (A14) 



ness coefficient k - Obviously, k is identified with lead us to attain the stiffness coefficients 



oi 0) \ r 
J 0) = 

u N — 

we obtain kq = —Em/[({ + vm)tn\- Moreover, the uni- 
form contraction energy UjfP should be represented as 



■ N associated with the unit radial displacement fen 
1. Hence, imposing an appropriate value of ao, 



Or \r=r N ,5p, N =X,S9 N =o/ COS 7l9 and fci2 



-o r / | r=rjVi(5/ZN . = o i 5p N= i/ cos n9 as 



U 



(0) 



(o)\- 

'JV J 



2(n + l)(l-^,)-l F 
fc n = /-, ttt; : ^r^M, 



r^d9. Consequently, we ob- 



tain the explicit form of U M that depends on Em , vm 



(1 + i/ju )(3 - 4j/ m ) 
2(n + 1)(1 - i/ M ) - n 
(1 + ^ M )(3 - 4^/) 



(A15) 
(A16) 



8 



and &22 = fell) &2i = ^12- Substituting the results (|A13[1 - 
(|A16[) into Eq. ([5]), we finally obtain the explicit form of 



AU 



(n) 
M • 
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